// global MTBQ score (NB. missing data must be coded as .)
* generate new variable where 0 = not difficult/does not apply, 1 = a little difficult, 2 = quite difficult, 3 = very difficult, 4 = extremely difficult
recode p1_ms_6_1_1difmeds (1 0 = 0) (1 = 0) (2 = 1) (3 = 2) (4 = 3) (5 = 4) , generate (p1_ms_6_1_1difmedsS)
recode p1_ms_6_1_2remem (1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_2rememS)
recode p1_ms_6_1_3cost (1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_3costS)
recode p1_ms_6_1_4collct(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_4collctS)
recode p1_ms_6_1_5monitr(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_5monitrS)
recode p1_ms_6_1_6arr_appt(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_6arr_apptS)
recode p1_ms_6_1_7dif_hcp(1 0 = 0)(1 = 0) (2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_7dif_hcpS)
recode p1_ms_6_1_8att_appt(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_8att_apptS)
recode p1_ms_6_1_9accss(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_9accssS)
recode p1_ms_6_1_10comm_hlp(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_10comm_hlpS)
recode p1_ms_6_1_11info (1 0 = 0)(1 = 0) (2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_11infoS)
recode p1_ms_6_1_12chang(1 0 = 0) (1 = 0)(2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_12changS)
recode p1_ms_6_1_13famhlp(1 0 = 0)(1 = 0) (2 = 1) (3 = 2) (4 = 3) (5 = 4), generate (p1_ms_6_1_13famhlpS)

* global MTBQ score if < 50% items (<6 questions) are missing we took the average of the scores, excluding items 3, 9 and 10
egen globalMTBQ=rmean( p1_ms_6_1_1difmedsS p1_ms_6_1_2rememS p1_ms_6_1_4collctS p1_ms_6_1_5monitrS p1_ms_6_1_6arr_apptS p1_ms_6_1_7dif_hcpS p1_ms_6_1_8att_apptS p1_ms_6_1_11infoS p1_ms_6_1_12changS p1_ms_6_1_13famhlpS )if missing(p1_ms_6_1_1difmedsS) + missing(p1_ms_6_1_2rememS) + missing( p1_ms_6_1_4collctS) + missing( p1_ms_6_1_5monitrS) + missing(p1_ms_6_1_6arr_apptS) + missing(p1_ms_6_1_7dif_hcpS) + missing(p1_ms_6_1_8att_apptS) + missing(p1_ms_6_1_11infoS) + missing(p1_ms_6_1_12changS) + missing(p1_ms_6_1_13famhlpS) < 6 

generate globalMTBQ1 = globalMTBQ*25

summarize globalMTBQ1
tabulate globalMTBQ1
histogram globalMTBQ1

* global MTBQ score if < 50% items (<7 questions) are missing we took the average of the scores, including items 3, 9 and 10
egen globalMTBQ=rmean( p1_ms_6_1_1difmedsS p1_ms_6_1_2rememS p1_ms_6_1_3costS p1_ms_6_1_4collctS p1_ms_6_1_5monitrS p1_ms_6_1_6arr_apptS p1_ms_6_1_7dif_hcpS p1_ms_6_1_8att_apptS p1_ms_6_1_9accssS p1_ms_6_1_10comm_hlpS p1_ms_6_1_11infoS p1_ms_6_1_12changS p1_ms_6_1_13famhlpS )if missing(p1_ms_6_1_1difmedsS) + missing(p1_ms_6_1_2rememS) + missing(p1_ms_6_1_3costSmissing) + missing( p1_ms_6_1_4collctS) + missing( p1_ms_6_1_5monitrS) + missing(p1_ms_6_1_6arr_apptS) + missing(p1_ms_6_1_7dif_hcpS) + missing(p1_ms_6_1_8att_apptS) + missing(p1_ms_6_1_9accssS) + missing(p1_ms_6_1_10comm_hlpS) + missing(p1_ms_6_1_11infoS) + missing(p1_ms_6_1_12changS) + missing(p1_ms_6_1_13famhlpS) < 7 

generate globalMTBQ1 = globalMTBQ*25

summarize globalMTBQ1
tabulate globalMTBQ1
histogram globalMTBQ1